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Abstract 

The kinematic properties of tidal debris from an orbiting Galactic satellite is presented, on the assump- 
tion that its central part once contained the most massive Galactic globular cluster, u Cen. We simulate 
dynamical evolution of a satellite galaxy that follows the present-day and likely past orbits of uj Cen and an- 
alyze the kinematic nature of their tidal debris and randomly generated Galactic stars comprising spheroidal 
halo and flat disk components. It is found that the debris stars show a retrograde rotation at ~ —100 km s~^, 
which may accord with a recently discovered stellar stream at radial velocity of ~ 300 km s~^ towards the 
Galactic longitude of ~ 270°. These stars also contribute, only in part, to a reported retrograde motion of 
the outer halo in the North Galactic Pole direction, without significantly modifying local halo kinematics 
near the Sun. The prospects for future debris searches and the implications for the early evolution of the 
Galaxy are briefly presented. 
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1 Introduction 

It is well known that the most massive globular cluster 
in the Milky Way, u Cen, shows unique properties in its 
metallicity content, internal kinematics, and structure. For 
instance, ui Cen shows a wide spread in metallicity un- 
like other Galactic globular clusters (e.g. Norris, Freeman, 
& Mighell 1996): its metallicity distribution is peaked at 
[Fe/H]~ —1.6, along with a second smaller peak at [Fe/H]~ 
— 1.2 and a long tail extending up to [Fe/H]~ —0.5. Also, 
the metal-rich stars in u Cen are largely enhanced in a- 
process elements relative to those in globular clusters and 
field stars with similar metallicities (e.g. Norris & Da Costa 
1995). This suggests that the ejecta from low-mass, asymp- 
totic giant branch (AGB) stars had to be retained and 
incorporated into the next-generation stars. However, al- 
though w Cen is most massive (5 x 10^ Mq), it is unable to 
retain the AGB ejecta, as shown by Gnedin et al. (2002). 
Thus, an isolated formation of u Cen is unlikely, because 
the enriched gas would easily be lost by encountering the 
Galactic disk. The most viable explanation for the unique- 
ness of oj Cen is that it was once the dense nucleus of a 
dwarf galaxy, i.e., a nucleated dwarf (Freeman 1993). A 
gravitational potential of progenitor's dark matter would 
help retaining the enriched gas and let the cluster being 
self-enriched at least over a few Gigayears. 

In this contribution, we pursue the possible kinemati- 
cal evidence for the existence of such a dwarf galaxy, which 
was already disrupted in the past. Dinescu (2002) first in- 
vestigated this issue, by examining the possible signature 
of the progenitor's tidal debris among nearby metal-poor 
stars in the catalog of Beers et al. (2000, BOO). She identi- 
fied a group of stars with —2.0 <[Fe/H]< —1.5, which de- 
parts from the characteristics of the inner Galactic halo but 
has retrograde orbits similar to uj Cen. Her simplified dis- 
ruption model of the progenitor galaxy demonstrated that 
trailing tidal debris, having orbital characteristics similar 
to the cluster, can be found in the solar neighborhood, al- 
though the concrete spatial distribution and kinematics of 



the debris stars remain yet unclear. 

Dinescu's work motivates us to undertake a more re- 
fined approach to the issue, i.e., to conduct an N-body sim- 
ulation for the tidal disruption of u) Cen's progenitor galaxy 
(Mizutani et al. 2003). We obtain the characteristic struc- 
ture and kinematics of its debris stars and compare with 
various observations showing signatures of recent merging 
events in the Milky Way (Gilmore, Wyse, & Norris 2002, 
GWN; Kinman et al. 2003, K03; Chiba & Beers 2000, CB). 
In particular, we show that a recently identified stream of 
stars at radial velocity of ~ 300 km s~^ (GWN) is a natural 
outcome of the current disruption model, without signifi- 
cantly modifying local halo kinematics near the Sun. 



2 Model of a Progenitor Galaixy 

For the purpose of enlightening the characteristics of tidal 
debris, we first assume a fixed external gravitational poten- 
tial (representing the Milky Way) for the dynamical evo- 
lution of an orbiting dwarf galaxy. The potential consists 
of a spherical Hernquist bulge ^^(r), a Miyamoto-Nagai 
disk <^d(R,z), and a logarithmic dark halo $h(r), where r 
is the Galactocentric distance and [R, z) denote cylindri- 
cal coordinates. Each is given as, ^^(r) = —GM^jjir -I- a), 
$d(i?,z) = -GMi/^R^ + {h + \jz^ + c2)2, and $h(r) = 
Wh/21n(r^ -h d^), where Mb = 3.4 x 10^" Mq, a = 0.7 kpc. 
Mi = 10" Mq, 6 = 6.5 kpc, c = 0.26 kpc, vn = 186 kms"S 
and d = 12 kpc. This choice yields a circular velocity of 
228 km s~^ at the solar circle of Rq = 8 kpc and a flat ro- 
tation curve outside Rq. On our adoption of a logarithmic 
dark halo, we note that while a dark halo may be modeled 
by the universal profile proposed by Navarro, Frenk, and 
White (NFW profile) , such a profile as obtained by N-body 
simulations can be modified by the inclusion of baryonic 
gas, mostly in the inner part of a galaxy where its density 
is increased owing to baryonic infall, in such a manner that 
the density profile may be changed to an isothermal one 
(e.g. Barnes 1987). It is also noted that the adoption of 
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the NFW profile does not affect our results, as the iruier 
density profile in concern is dominated by the bulge and 
disk components. A dwarf galaxy orbiting in this galactic 
potential is then represented by self-gravitating particles 
following a King model, where the central density, cen- 
tral velocity dispersion, and core radius are given as 0.3 
Mq pc~^, 18.1 km s~^, and 0.56 kpc, respectively. In ad- 
dition, a particle with the mass of 5 x 10^ Mq representing 
uj Cen is placed at the center of the galaxy. This setting 
yields the total mass of the system as Mtot = 5.79 x 10* 
Mq. On the other hand, the total mass of stars in this 
hypothetical galaxy is estimated from the mean metallicity 
of stars in uj Cen (([Fe/H]) ~ —1-6), combined with the 
metallicity-luminosity relation for the Local Group dwarfs 
(Cote et al. 2000) and the mass-to-light ratio (assuming 
M/L ~ 4), yielding M^tars ~ lO'' M©. Thus, the mass of 
the simulated galaxy is largely dominated by dark matter. 
The model galaxy is represented by a collection of 10* parti- 
cles and the self-gravity is calculated in terms of a multiple 
expansion of the internal potential to fourth order (White 
1983; Zaritsky & White 1988). 

This model dwarf is disrupted by Galactic tides in the 
course of its orbital motion, whereas its dense core is ex- 
pected to survive and follow lj Ceri's orbit. We examine 
two representative orbits for the progenitor, models 1 and 
2: model 1 follows the current orbit of lo Cen, whereas for 
model 2, we calculate an orbit back to the past over ~ 2 
Gyr from its current position and velocity taking into ac- 
count dynamical friction and then set a progenitor galaxy 
on its non-decaying orbit. These two models provide us 
with satisfactory information on the generic properties of 
a tidally disrupted progenitor and we postulate that the 
realistic nature of their debris is midway between these 
model predictions. We calculate u! Cen's orbit, based on 
the current distance D — 5.3 ± 0.5 kpc from the Sun, po- 
sition {l,b) = (309°, 15°), proper motion cos S,ns) = 
(-5.08 ± 0.35,-3.57 ± 0.34) mas yr"\ and heliocentric 
radial velocity vios = 232.5 ± 0.7 km s^^ (Dinescu, Gi- 
rard, & van Altena 1999, DGvA). This orbit for model 1 
is characterized by frequent disk crossings with a period of 
Torb = 0.8 X 10* yr, retrograde motion, and apo and peri- 
centric distances {rapo,rperi) ~ (6.4, 1.1) kpc. For model 2, 
we obtain Torb ~ 1.5 x 10* yr and {rapo,rperi) ~ (11.3, 3.0) 
kpc. In both experiments, we place a progenitor galaxy at 
apocenter to maximize its survival chances. 

We plot the spatial distribution of the tidally disrupted 
debris in Figure 1. Upper (middle) panel shows model 1 
(model 2) after the 1.37 (1.86) Gyr orbital excursion of 
the progenitor gala:xy. Lower panel shows the orbit of the 
galaxy center. A rosette-like feature of the debris becomes 
steady after about eight orbital periods. Model 1 results in 
more compact distribution than model 2, which reflects the 
difference in orbital radii. Figure 2 shows the velocity dis- 
tributions of the debris particles in cylindrical coordinates 
{vr,v^,Vz). As is evident, both models provide essentially 
the same debris kinematics: most notorious is a sharply 
peaked distribution at ~ —100 km s^^, arising from a 
retrograde orbit of a progenitor. These kinematics suggest 
that the difference in model 1 and 2 resides only in the 
spatial extent of the debris. 
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Figure 1: Upper (middle) panel shows the spatial dis- 
tribution of the tidally disrupted debris for model 1 
(model 2) after the 1.37 (1.86) Gyr orbital excursion 
of to Cen's progenitor galaxy. Lower panel shows the 
orbit of the galaxy center for model 1 (solid line) and 
2 (dotted line). The plots are projected onto three or- 
thogonal planes, where the Sun is located at x = — 8 
kpc and xy corresponds to the disk plane. A frame 
measures 15 kpc on a side of each panel. The current 
position of w Cen is at {x, y, z) = (—4.8, —4.0, 1.4) kpc. 



3 Kinematics of Stream and Galactic Stcirs 
3.1 Method 

We analyze the kinematics of both the simulated debris and 
other Galactic stars generated randomly by a Monte Carlo 
method. The metal-poor halo is modeled as a flattened 
spheroid p oc (i?^ + z"^ /q^)^'^'"'^^ , where q is an axis ratio 
ranging 0.55-0.7, anisotropic velocity ellipsoid (ur, a^, Oz) ~ 
(154,121,96) km s~^, and small mean rotation (v^) = 24 
km s^^, as found for halo stars with [Fe/H]< —2 near 
the Sun (CB). Thin and thick disks are modeled as p oc 
exp{—R/Rd)sec'^{z/zd), where Ra = 3.5 kpc and za — 0.3 
(1) kpc for thin (thick) disk. Both disks rotate at 200 
km s~^, having velocity ellipsoids of (34,25,20) km s~^ 
and (46, 50, 35) km s~^ for thin and thick disks, respectively 
(CB). The relative fraction of each component is fixed using 
observed local densities near the Sun, in such a manner that 
the halo and thick-disk densities at -D < 1 kpc are 0.2 % 
and 2 % of the thin-disk density, respectively (Yamagata & 
Yoshii 1992). 

In our model of lj Cen's progenitor galaxy, the self- 
gravitating particles represent both stars and dark matter 
particles. We note that the correct conversion of the mass of 
the simulated system into that of the presumed stellar sys- 
tem alone remains yet uncertain, because the correct M/L 
ratio of stars as well as the amount of dark matter in the 
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Figure 2: Velocity distribution of the debris particles 

in cylindrical coordinates, vr (dotted), v,p (solid), Vz 
(dashed), for model 1 (a) and 2 (b). 



progenitor galaxy is unavailable. As a useful model param- 
eter to incorporate this ambiguity for the current kinematic 
analysis, we set a quantity / as the fraction of the debris 
particles relative to halo stars near the Sun, so that the nor- 
malization of the halo density is fixed for the given number 
of the neighbor debris and /. 

A typical value of / for the conversion of the simulated 
particles to the stars is estimated in the following manner. 
Model 1 (model 2) yields 21 (74) particles at D < 2 kpc, giv- 
ing the mass density of Pg = 0.4(1.3) x 10"'' Mq pc"^ near 
the Sun, whereas the total mass density and metal-poor 
halo density have been derived as 8x 10~^ M© pc~^ (Gates, 
Gyuk, & Turner 1995) and 6.4 x 10"^ M© pc"^ (Gould, 
Flynn, & Bahcall 1998), respectively. Then, if the debris 
stars (with M stars ~ 10^ Mq) are distributed in the same 
manner as the simulated particles (with Mtot ~ 5.74 x 10* 
Mq), which would be a reasonable approximation in view 
of the dissipationless nature of stars, the mass density of 
the debris stars in the solar neighborhood can be estimated 
as (Mstara/Mtot)pg = 0(10"*^) Mq pc"^, which is about 
1 % of the halo density. Tims, /', defined here at f < 2 
kpc, is expected to be of order of a few percents. 



3.2 Radial velocity distribution at I ~ 270° 

GWN recently reported their spectroscopic survey of ~ 
2000 F/G stars down to V = 19.5 mag, in the direction 
against Galactic rotation {I, h) = (270°, -45°) and (270°, 4-33 
where radial velocities in combination of distances largely 
reflect orbital angular momentum. Their calibration of the 
stars at distances of a few kpc from the Sun leads to the 
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Figure 3: (a) Distribution of the heliocentric radial 
velocities in the direction of the GWN survey, for the 
debris stars of model 2 and randomly generated halo 
stars with q — 0.7. We select the stars at 1 < I? < 
5 kpc and 260° < I < 280° in the fields of 30° < 
|6| < 50° (with / = 1 % and 3 % for dotted and solid 
histograms, respectively) and 20° < |6| < 40° (with 
/ = 3 %: dashed histogram), (b) The same as (a) but 
incorporating the randomly generated disk stars as well 
for 30° < |6| < 50° with / = 5 % (solid histogram). 
Dotted histogram denotes the contribution from the 
metal- poor halo alone. 



discovery of two stellar streams at vios ~ 100 km s~^ and 
~ 300 km , which are not explained by any existing 
Galactic components (see their Fig. 2 and 3). While the 
stream at vios 100 km ~^ was reproduced by their model 
of a merging satellite in prograde rotation, the stream at 
Vios ~ 300 km ~^ remains yet unexplained. 

Figure 3a shows the vios distribution for the simulated 
debris of model 2 and halo stars oi q = 0.7 (i.e. without 
disks) at 1 < D < 5 kpc, 260° < / < 280°, and two fields 
for 6. Figure 3b shows when disk stars are incorporated. As 
is evident, the debris stars from w Cen's progenitor galaxy 
form a local peak at vios ~ 300 km s~^, which is provided 
by many stars having ~ —100 km s^^ (Fig. 2). This 
is in good agreement with the vios ~ 300 km s~^ stream 
discovered by GWN. A more flattened halo than the case 
q — 0.7, which is possible within a range of 0.55 < q < 0.7 
(CB), yields a higher peak, since the density contrast of the 
debris relative to halo is made higher in this survey region. 
°llVe also make a plot for model 1 but do not present it here 
as it is essentially the same as Figure 3 for model 2. This 
reflects that model 1 shows the same velocity distribution as 
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model 2 (Fig. 2), although to attain the same peak height 
at ~ 300 km s~^, / be a few factor larger and the selected 
range of / be a few degree higher than the respective values 
in model 2, because of less number of debris particles near 
the Sun. This rule applies to other considerations below as 
well. 



3.3 Kinematics at the North Galactic Pole 

The halo kinematics at the North Galactic Pole (NGP) de- 
serve special attention. Majewski (1992) suggested that 
the outer halo at the NGP shows a retrograde rotation 
(vtf,) ~ —55 km s^^ at 2: > 4 kpc. Also, K03 reported 
from their analysis of RR Lyrac and blue horizontal branch 
stars that the halo at 2 < z < 12 kpc shows a retrograde 
rotation at (v^,) ~ —65 km s~^. On the other hand, halo 
stars near the Sun show no significant retrograde rotation 
(CB). 

For the purpose of elucidating the effect of the debris 
stars on the NGP kinematics, we select the debris of model 
2 and randomly generated stars at 6 > 70° and 2 < D < 5 
kpc (resembling KOS's selection). Since the observational 
determination of full space velocities involves rather inaccu- 
rate information of proper motions compared to radial ve- 
locities, we convolve the velocity distribution of stars with 
a Gaussian distribution for velocity errors, having 1 a of a 
typical 30 km s^^ error. Figure 4 shows the velocity dis- 
tributions at various values of /. It follows that the debris 
stars provide a non-Gaussian feature in the velocity distri- 
bution. The distribution shows an extra peak at ~ —100 
km s~^, whereby the debris stars decrease {v^) from that 
of a pure halo sample. However, the change of {v^} by the 
inclusion of the debris stars with / = 5 % amounts to only 
-19 (-14) km s"^ for q = 0.55 (0.7), which are still in- 
sufficient for explaining the reported (v^) = —35 ~ —75 
km s^"*^. For the vr and Vz distributions, their means at 
the cases of q = 0.7 and / = 5 % are slightly increased to 
(vr) = 7 km s~^ and (vz) = 3 km s~^, respectively, whereas 
the velocity dispersions are decreased by the amounts of 6 
km s~^ and 5 km s~^, respectively, when the debris stars 
are considered. Note that if we extend the selection of the 
stars at higher z or instead consider model 1, the changes 
of (v^) become smaller than the above mentioned values, 
because there are no debris stars in our current model (Fig. 
1). Thus, it is safe to conclude that the debris stars con- 
tribute only in part to a reported retrograde motion at the 
NGP. 



3.4 Local halo kinematics 

We select the simulated stars at the cylindrical coordinates 

of 6.5 < R < 9.5 kpc and z < A kpc and at the distance from 
the Sun of _D < 4 kpc (as was drawn by CB), convolve the 
velocities with a Gaussian error distribution of 1 cr = 30 
km s~^, and compare with the corresponding staxs with 
[Fe/H]< —2 in BOO. It follows that the non-Gaussian fea,- 
ture of the velocity distribution is much weaker than that at 
the NGP: the change of (w^) for / = 5 % is only -9 (-10) 
km s~^ in model 1 (model 2). This is due to the charac- 
teristic debris distribution, where the number density near 
^; = is smaller than at high \z\ (Fig. 1). 
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Figure 4: Velocity distribution of the simulated stars 
of model 2 in the North Galactic Pole direction. In this 
plot, wc select the stars at & > 70° and 2 < £> < 5 kpc 
and convolve their derived velocity distribution with a 
Gaussian error distribution having 1 u of 30 km s~^. 
The axis ratio of the mctal-poor halo is set as g = 0.7. 
Dotted, solid, and dashed lines denote / = 1 %, 3 %, 
and 5 %, respectively. 



4 Discussion 

Our simple model of an orbiting dwarf galaxy that once 
contained a» Cen predicts a sequence of tidal streams in 
retrograde rotation and their existence can be imprinted 
in kinematics of nearby stars, especially in the direction 
against Galactic rotation (GWN) and at the NGP (K03), 
while local halo kinematics remain unchanged. The sim- 
ulated streams are mostly distributed inside the solar cir- 
cle, as suggested from the current orbital motion of uj Cen 
(DGvA; Dinescu 2002). In contrast to Sgr dwarf galaxy 
having polar orbit (Ibata ct al. 1997), the orbit of Cen's 
progenitor galaxy is largely affected by a non-spherical disk 
potential, where the orbital plane exhibits precession with 
respect to the Galactic Pole, causing self-crossing of tidal 
streams in the disk region (Fig. 1). The projection of the 
orbit perpendicular to the disk plane shows an 'X'-likc fea- 
ture, thereby leaving denser streams at high \z\ than at low 
\z\ for a given radius. These characteristic spatial distribu- 
tions of the debris stars give rise to more significant effects 
of the debris at the NGP than in the solar neighborhood, 
as shown here, although the current simulation failed to 
reproduce the reported largely retrograde rotation at NGP 
from the uj Cen debris alone; perhaps, other, yet unknown 
halo substructures must be considered to reproduce the ob- 
servations. 

Existing kinematic studies of Galactic stars to search 
for a signature of uj Cen's progenitor galaxy are yet con- 
fined to nearby stars, where the significance of the debris 
streams is modest, as shown here. Searches of stars inside 
the solar circle are more encouraging (Fig. 1), in particular 
in the directions of I ~ 320° and I ~ 50° , where we expect 
the presence of high-velocity streams at vios = 200 ~ 300 
km s~^ and —400 ~ —300 km s~^, respectively. Future ra- 
dial velocity surveys of these fields including the sample of 
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the Sloan Digital Sky Survey or planned Radial Velocity Ex- 
periment are worth exploring in this respect. Also, detailed 
abundance studies of candidate stream stars will be intrigu- 
ing, because such stars may exhibit different abundance 
patterns from field halo stars, as found in dwarf galaxies 
(Shetronc, Goto, & Sargent 2001). 

In this work, we adopt a fixed external gravitational 
potential for the calculation of an orbiting dwarf galaxy 
and its dynamical evolution, thereby neglecting dynami- 
cal friction against the satellite. However, to set more re- 
fined limits on its dynamical history, it is required to fully 
take into account frictional effect on the orbit as well as 
dynamical feedback of a satellite on the structure of the 
Galactic disk. For instance, Tsucliiya, Dincscu, and Kor- 
chagin (2003) reported their numerical models for a satel- 
lite with strong orbital decay and succeeded to reproduce 
the current orbit of u} Gen from the launch of the pro- 
genitor satellite at 58 kpc from the Galactic center. They 
also showed that the debris particles at 3 Gyr are already 
smeared out and distributed in a flattened disk without 
having significant stream-like features as obtained here. It 
is noted that whether or not stream-like features survive 
by the current epoch depends on when and how a satellite 
galaxy is disrupted by Galactic tides, or in other words, the 
observational information on such features and the compar- 
ison with simulation results are useful for placing impor- 
tant constraints on when a satellite merging occurs. For 
this purpose, we are currently undergoing full N-body sim- 
ulations of several host-satellite systems using GRAPES, 
where both host and satellite galaxies are represented by 
lively dark halos and stellar components. Our goals with 
the use of GRAPES are to set tight limits on the initial to- 
tal mass, internal mass distribution, and orbital motion of a 
progenitor galaxy, as well as the timing of merging with the 
Galactic disk, based on the comparison between the obser- 
vation (e.g., RAVE) and simulation results for stream-like 
structures in the Milky Way. More details will be reported 
elsewhere (Mizutani & Ghiba, in preparation). 
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